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Abstract 

We investigate the thermodynamical content of a cellular model which describes 
nuclear fragmentation as a process taking place in an excited disordered system. The 
model which reproduces very well the size distribution of fragments does not show 
the existence of a first order phase transition. 
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1. Introduction 

The description and understanding of nuclear fragmentation has remarkably 
progressed during these last years. There remains nevertheless a long standing 
question for which no completely satisfactory answer has yet been found. It concerns 
the existence of an energy triggered signal which would reveal the existence of a 
phase transition in nuclear matter. 

Among the advances which have come out through the confrontation of experi- 
mental measurements [1-5] with theoretical investigations one may mention the 
behaviour of final fragment charge distributions. They show that the fragmentation 
process leads to universal properties which can be interpreted in the framework of 
very simple descriptions [6] . In particular, the interpretation of percolation analyses 
points towards the existence of a second order phase transition which takes place in 
an excited disordered system [7-8]. More recently, experimental data [9] were used 
in order to relate the excitation energy to a temperature which characterises the 
fragmenting system produced in peripheral collisions. The caloric curve which comes 
out of these data shows a typical plateau in temperature, with some kind of up-swing 
on the lower side of the energy interval over which the effect is observed. Whether 
this result and its interpretation will be definitely confirmed or not is, at present, 
an open question. First, it has to be cross-checked that the temperature which was 
obtained by means of a model [10] is consistent with other experimental methods 
aimed to fix this quantity. Up to now, this is not clear [11]. Second, the concept of 
temperature itself makes sense only if the system is in thermodynamical equilibrium 
which is tacitly supposed to be realised, although it can never be rigorously so. 
Indeed, one expects the presence of more or less strong non-equilibrium features 
like collective flow, even in peripheral kinematics. This point has to be kept in 
mind when a comparison between theoretical model calculations and experimental 
results is performed. 

Many different theoretical models have been proposed in order to describe 
nuclear fragmentation. The specific question concerning the thermodynamic pro- 
perties of fragmenting nuclei and the possible existence of different phases has been 
investigated by means of time-dependent approaches [12-14] as well as thermody- 
namic equilibrium models which describe the freeze-out stage of the process in the 
framework of the canonical [16] and the micro canonical ensembles [15]. Both of 
these last calculations lead to a caloric curve which shows the characteristic be- 
haviour of a first order phase transition taking place in a finite system. This may 
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be interpreted as a support of the experimental result, although, as already stated, 
it is not clearly established that the physical conditions which are experimentally 
realised are those introduced by the models. 

In a former work [7,8] we proposed a classical microscopic description in which 
a fragmenting nucleus is an excited and disordered system. The approach is able 
to reproduce the results of bond percolation as well as the ALADIN data [17,18] 
concerning charge distributions of fragments and related observables. This success 
legitimates the interest in further phase space properties of the model. In particular, 
it is worthwhile to see whether it can or cannot reproduce the thermodynamic 
properties predicted by the recent experiments. 

The organisation of the present paper is the following. In section 2 we first recall 
the essential features of our approach and describe the way we construct a partition 
function of the system in the framework of the microcanonical ensemble (fixed 
energy) . We discuss the way through which we establish the correspondence between 
the energy and a temperature and present the results of our numerical investigations. 
In section 3 we develop an analytical version of our approach which we confront with 
the numerical findings. We discuss our results and draw conclusions in section 4. 

2. Description of a disordered and excited system of particles and its 
thermodynamical properties 

2.1 Cellular model approach to nuclear fragmentation 

We start with the description of the model which has already been introduced 
in Refs. [7,8]. We consider a finite compact volume V in 3-dimensional space, in 
practice a cube. This volume is divided into A equal cubic cells of linear dimension 
d. Each cell is occupied by a classical particle, i.e. V = Ad 3 . Each particle i is 
characterized by its phase space coordinates, i.e. its position f*i with respect to 
the centre of mass of the system and its linear momentum p^. The coordinates 
(xi, yi, Zi) are chosen as Xj = Xi +ed/2 and similarly for j/j and Zi. The coordinate 
fi (xi , yi , Zi ) corresponds to the centre of the cell and £ is a random number drawn 
from a distribution which, in practice, is chosen to be uniform in the interval [0,1]. 
The momenta Pi are drawn randomly and uniformly in a sphere of arbitrary large 
radius p max - We introduce the maximum kinetic energy of a particle as E max = 

The particles interact by means of two-body interactions. In practice we in- 
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troduce a strong short range nuclear potential. For reasons of consistency with 
former calculations [8] we essentially use the interaction of Ref. [19] which does 
not explicitly differentiate protons and neutrons. If we select another potential in 
the sequel we shall mention it in due place. We choose part of the nucleons to be 
protons which we distribute uniformly in cells over the whole space and implement 
the Coulomb potential which acts between charged particles. The potential energy 
of the system is obtained by adding up the interaction energy between the parti- 
cles. We have also considered the case where the nuclear interaction is restricted 
to nearest neighbours. This approximation changes absolute values of energies but 
does not affect the conclusions which will be drawn below. 

Former use of this classical approach has nicely proven that the present des- 
cription of excited fragmenting nuclear systems is able to reproduce the fragment 
size distributions and related observables which are obtained in the framework of 
ordinary bond percolation [7,8]. The calculations show that the results are very 
robust with respect to the explicit choice of the interaction, the shape of the external 
surface of the system and even the criteria which are used in order to define bound 
systems of particles, i.e. fragments. 

2.2 Thermodynamical aspects 

2.2.1 Motivations 

The agreement between percolation calculations which are generic for the des- 
cription of disordered systems and experimental results suggests a possible inter- 
pretation of fragmentation as a critical phenomenon, more precisely a second order 
phase transition in presence of finite size effects since it takes place in a small sys- 
tem. A recent interpretation of experimental data raises a new question. Indeed, 
the determination of the temperature of the excited fragmenting system and the 
knowledge of the excitation energy of projectile-like nuclei formed in peripheral col- 
lisions lead to a caloric curve which exhibits the characteristic plateau corresponding 
to latent heat generation. This effect which is seen in infinite systems undergoing 
a first order phase transition is interpreted as a liquid-gas-like transition [9]. This 
result which as already said, is still under critical investigation and its interpretation 
led us naturally to the question as to whether our approach can reproduce such an 
effect or not. This is the point we want to present and discuss in the sequel. 

2.2.2 Generation of microcanonical configurations 

The total number of configurations which can be generated for a fixed energy 
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E, number of particles A and volume V = Ad 3 is designated by Z (E, A,d), the 
micro canonical partition function of the system when it is in thermo dynamical 
equilibrium. 

In practice we construct this quantity numerically in the following way. We 
generate N configurations with our cellular model as described above and calcu- 
late the potential energy of the system by summing up the two-body interaction 
energy which acts between all pairs of particles of the system, i.e. the nuclear and 
the Coulomb contributions. The number Z of charged particles is such that Z/A 
corresponds to nuclei energetically located in the valley of stability. Doing this way 

(k) 

we obtain N configurations with potential energy {E pot , k = 1, . . . , N}. 

In a second step we generate a set of configurations for which we determine 
the kinetic energy content. For each configuration we choose randomly p max (see 
above) out of a uniform distribution in an interval [0, P up ] where P up is an upper 
limit for the momenta which we fix arbitrarily. We then draw randomly, for fixed 
p max , an ensemble of momenta = 1,...,A} in the Fermi sphere of radius 

Pmax as already described above. Proceeding this way we generate a set of events 
with total kinetic energy {E^ n , k = 1, . . . , N}. The number of events is generally 
of the order of 10 6 . 

We fix then the total energy of the system to lie in bins [E, E + AE] and 
count the number of configurations which correspond to a total energy E = Ek in + 
E P ot ly m g hi this bin by means of the combination of the kinetic and potential 
energy configurations which were generated as described above. The number of 
configurations which are found gives the partition function Z (E, A,d). This is 
done for all energies which the system can reach, from a minimum value E m i n to 
the maximum which is E up = A ■ P^ p /2m. E up is necessarily a finite quantity 
because it is necessary to define a finite upper bound in numerical simulations. 
This leads to the fact that Z (E, A, d) may start to decrease for the higher values 
of E. This is an unphysical effect. Hence we must restrict the further use of the 
partition function to an interval of energy which does not take its high energy tail 
into account. 

2.2.3 Thermodynamical analysis of the disordered system. Results 

The knowledge of the microcanonical partition function leads to the entropy 
S (E, A, d) and the temperature through 



T 



-l 



(ds/dE) v 
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which establishes the relationship between T and E, i.e. the so called caloric curve. 

We have determined Z (E, A, d) for different values of d, i.e. different total 
volumes. In order to work out the derivative of S = In Z we have introduced an 
analytical parametrisation of S over the whole range of values in terms of Chebyshev 
polynomials of different orders. Using statistically significant samples of configura- 
tions (of the order of 10 6 ) it is possible to approximate S with an excellent precision 
(x 2 < !)• Fig. la shows the behaviour of Z obtained by means of the numerical 
simulations described above and the result of the analytical fit which is practi- 
cally indistinguishable from the numerical result. The corresponding caloric curve 
is represented in Fig. lb. There is no sign for the existence of a plateau in the 
temperature which would signal a first order phase transition. 

In Fig. 2 we present the behaviour of (dS/dE)^} 1 for different values of V 
which correspond to values of d in the interval [1.8, 4] fm (average densities p in 
the interval [0.016, 0.17] fm -3 ). One observes a strong dispersion of the curves as a 
function of the energy which indicates that this quantity is sensitive to the volume 
in which the system moves. The large negative energy thresholds come from the 
fact that the nucleon potential we chose [19] and which was also used in Refs. 
[7,8] overbinds heavy systems (here A = 216). In our former study of fragment size 
distributions [7,8] we mixed up events which correspond to systems of different sizes. 
Since, for different experimental events, the nucleus may break up into fragments at 
different densities it is tempting to collect in a uniform way events corresponding to 
different values of d. This is what we did in order to construct an "entropy" which 
is averaged over different values of d, < S >d- In Fig. 2 we show the behaviour 
of (d < S >d jdE)y X (dotted line). Again, this quantity shows a monotonous rise 
with the energy. 

We also used the potential defined in Ref. [20] which differentiates protons and 
neutrons and hence leads to a quite different energy dependence of the system. As 
it can be seen in Fig. 3 one finds the same behaviour of the observables as in the 
former case. 

2.2.4 Comments 

We comment and discuss now the results obtained above. First, the fact that 
we take together events corresponding to different volumes violates (1) since dS/dE 
must be calculated for fixed volume. The only way out is to say that the quantity 
which corresponds to (d < S >a /dE)~ x is some effective temperature. It is 
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not clear whether it could be put in connection with temperatures related to the 
experimental situation where the fragmentation may go on in a variable volume. 

Second, it is interesting to confront the behaviour of the partition sum Z (E, A, d) 
with the multiplicity as a function of increasing energy. This is shown in Fig. 4. As 
expected there exists a clear correspondence between the increase of the number of 
configurations and the increase of the number of fragments and particles over the 
whole energy range. However this connection is not related to a phase transition, 
as already stated above (no exponential increase of Z (E, A, d) over a finite energy 
interval in which one would observe a coexistence of fragments and particles). 

Last, ther mo dynamical equilibrium which is a prerequisite for the use of (1). 
There is no reason why all the configurations we generate should correspond to such 
a situation since we do not implement a Metropolis Monte Carlo type algorithm in 
order to select those configurations which correspond to equilibrium. We try now 
to eliminate the problem concerning thermodynamical equilibrium by implementing 
an analytical approach. 

3. Analytic approach to the description of an excited disordered system 

As previously discussed, the numerical simulations which were performed above 
do not necessarily generate an equilibrated system. It raises the questions whether 
the outcome of our model would be different if this property would be satisfied. This 
should be done by implementing some Metropolis Monte Carlo algorithm. Such cal- 
culations may be done in the future. Here we try to work out an analytical model 
which guarantees thermodynamic equilibrium. Furthermore, its geometric and dy- 
namical properties are very close to those of the system which has been investigated 
above. We want to check the effects of possible deviations from equilibrium which 
is not realised in the calculations of section 2. 

3. 1 Derivation of the canonical partition function 

We construct the canonical partition function of a set of A classical particles 
in thermodynamic equilibrium at a temperature T = (3~ . The particles are dis- 
tributed in A identical cells of linear dimension d. As before, there is one particle 
in each cell. Then 

Z{P,A,d)= Yl e~ m{AA) (2) 

C({fi,Pi}) 
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TT f dpi f drie- pH ^ pi 
i= i Jv Pi Jv r . 



}) 



with 



V Pi and V ri are the volumes occupied by a particle in a phase space cell. We 
consider first a 1-dimensional system for which cells are segments of length d. The 
potential V is defined in terms of a truncated harmonic oscillator which acts between 
neighbouring particles % and j located at Xi, Xj such that 

V = oo \xi — xj\ < r c 

= Vq + C (xi - xj) 2 r c < \xi - xj\ < r 

= \Xi — Xj\ > To 

The parameter r c introduces a hard core. The effective range of the interaction 
r is chosen to be of the order of 2d. The strength Vq is negative and C = —Vq/tq. 

The potential V is aimed to simulate a realistic short range nuclear potential 
which is attractive and strongly repulsive at short distances. 

The model can be easily extended to 3 dimensions (see below). It is expected to 
be qualitatively close to the more realistic potentials used in numerical simulations. 
The long range Coulomb interaction is not taken into account. From our numerical 
experience this should not alter our conclusions. 

3.2 Kinetic and potential contributions to the partition function 

The kinetic contribution is trivially given by 

Z kin (p,A) = Y[ / d Pl e W (-(3p 2 /2m) = (27vm/(3) A / 2 (3) 
i=i J -°° 

The potential contribution can also be worked out explicitly. The potential 
energy reads 

V = AV + C [On - x 2 ) 2 + ...+ - x A ) 2 + (x A - Xl ) 2 } 

We impose periodic boundary conditions 

XA+l = Xi 
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and introduce new coordinates 



Xi = x® + 5xi ; Sxa+i = Sxi 

where X® defines the centre of the interval i. 
Then 

A A 

V = AV Q + ACd 2 + C ( 6x *+i ~ 5x i) 2 + 2Cd ( 6xi+1 ~ 6x ^ 

i=i i=i 

The periodic boundary condition eliminates the last term on the right hand 
side. Hence the potential contribution to the partition function reads 



Z pot (/?, A, d) = exp [-/3A(V Q + Cd 2 )]- J] / dSx k exp 

k=i J ~ b 



A rb 

^^(fei+l-fei)' 
i=l 



where [—6, +b] with b = d/2—r c is the range over which dxi can vary in the interval i. 

The exponent of the integrant in Z pot is a quadratic form which corresponds to 
a matrix with A equal diagonal elements 2/3C, two symmetrical subdiagonals with 
equal elements — f3C and two elements — (3C, one in the upper right corner and one 
in the lower left corner. This "tight-binding" matrix can be explicitly diagonalised 
and the eigenvalues are given by 



Xi = 2(3C (l - cos {^^j ^ = l3C\i 



(i = l,...,A) 



The lowest eigenvalue corresponds to = 0. The other eigenvalues are degen- 
erate of order 2, Xa-i = Xi. If A is even, the eigenvalue Xa/2 = 4/3C is the largest 
and it is non degenerate. 

It is also possible to work out the expression of the eigenvectors and to obtain 
the new boundaries in the new variables {5Xi} corresponding to the rotated system 
and associated with the eigenvalues {Xi}. 

For A odd one gets 

—b VA < 5Xi < + b VA corresponding to 

and 

— bi < 5X 2 i < + bi 

- Ci < SX 2 i+i < + Ci 
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where 



k=i 



2tt/c . 
cos I — — • I 



A-l 
2 



d = b 



k=i 



2tt/c 
sin I — — • % 



which correspond to the degenerate eigenvalue \. 
For A even the boundaries are 

-b VA < SX h2 < + bVA 

corresponding to Xa and Xa/2 
and 

— bi < 5X 2 i+i < + bi 

- a < 5X 2 i+2 < + c% 

corresponding to the degenerate eigenvalues Aj. Here bi and Cj are the same as 
those given above, except that the upper limit of the sums corresponds now to 
k = A/2 - 1. 

As a consequence of these developments the expression of Z po t can be explicitly 
integrated for finite A. For A odd one gets 

Z$(P, A, d) = exp [-pA{V + Cd 2 )] ■ 2b A 1 / 2 • (tt//3C) (a - 1)/2 



'(A-l)/2 

n 

i=i 



(A-l)/2 



II { erf [(PCX^b^ ■ erf [OSCA,) 172 ^] } (4a 



i=l 



and for A even the expression reads 

Z$ t (P,A,d) = exp [-(3A(V + Cd 2 )} ■ 2b A 1 / 2 • (tt/^CA^) 1 / 2 • (tv/(3C)^ 



'{A-2)/2 

n * 

i=i 



A/2-1 



J] {erf [(/^CA,) 1 / 2 ^] • erf [(^CA,) 1 / 2 ^] } (4b) 
- l 



where erf (x) is the error function and Aj = \/(3C (i = 1, . . . , A). 
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3.3 Approximate expressions 

It is now possible to write down Z (/?, A, d) as the product of Zki n (f3, A) given 
by (3) and Z pot ((3,A,d) given by (4a) or (4b). This is however of no practical 
use since it would not allow us to get an explicit expression of the microcanonical 
partition function Z (E 1 , A, d) which we would like to confront with our simulation 
results. 

For this reason we introduce an approximation in our calculation of Z po t. In 
practice, we replace the rotated hypercube which defines the volume in which the 
A — 1 rotated coordinates {SXi} vary by a hypersphere whose radius R is defined in 
such a way that its volume is the same as that of the original hypercube of dimension 
A — 1 (the coordinate corresponding to A^ = is not taken into account) . 



This procedure leads to 

R= [V ■T(A + l)/2] 1 / ( - A - 1) /n 1 / 2 

where V = (26) ^"^ • I ]7 A J for A odd. 

If we introduce this approximation we get the partition function 

z£> t (P, A, d) = exp [-pA(V + Cd 2 )] 2b A 1 ' 2 ■ (tt / pC)^/ 2 



(5) 



\A-l)/2 

n 



i=i 



(A-3)/2 



1 - exp(-/3CR 2 ) J2 



{(3CR 2 ) 



2\m 



m=0 



mi 



(6a) 



for the case where A is odd, and 



Z^ t (p, A, d) = exp [-PA(V + Cd 2 )] ■ 2b A 1 / 2 ■ (n/ pCX A/2 ) 1 / 2 ■ (tt/(1C) 



(A-2) 



■ erf 



bA^iPCXA/z) 1 



/2 



(A-2)/2 

n a . 



-1 r 



(A-4)/2 

1 - exp(-(3CR 2 ) 

m=0 



(f3CR 2 ) 



2\m 



ml 



(6b) 



for the case where A is even. In (6b) R is defined by (5), but A — 1 is replaced by 
A-2. 

The exact and approximate expressions Z pot and Z po t are close and it is easy 
to guess that their behaviour as a function of (3 is qualitatively similar. We have 
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checked this point by means of explicit numerical calculations in which we compare 
the product of error functions appearing in (4a) and (4b) with the R— dependent 
terms in (6a) and (6b). Calculations done for different values of A and different 
choices of parameters are shown in Fig. 5 for the quantities by which and Z (°) 
given by (4a) and (6a) differ (last terms of the expressions). The results show that 
these quantities have the same monotonous behaviour, with the same limit as (3 
goes to infinity. In both cases the derivatives of the curves go to zero as (3 goes to 
zero. 

3-4 Microcanonical partition function 

We use Z in order to construct the microcanonical partition function by means 
of an inverse Laplace transform 

Z (E, A, d) = — / d(3e? E Z(0, A, d) (7) 

Because of the explicit dependence of Z on (3 it is possible to work out Z explicitly. 

It should be noticed here that inverse Laplace transforms can be highly unstable 
quantities [21,22]. This fact is particularly troublesome when Z{(3) is numerically 
determined since any deviation of Z((3) can produce important deviations of Z(E) 
from its real value. In our case Z{(3) is known analytically. We believe that, even 
though Z(E) we obtain from this analytic expression is not the same as the exact 
one, Z(E) reflects the qualitative behaviour of the exact partition function Z{E) 
because of the analytical behaviour of both functions Z(f3) and Z((3) (see above 
and Fig. 5). This is sufficient in order to conclude confidently about the qualitative 
behaviour of E vs. T (see below and Figs. 6 and 7). 

Using (6a) it is straightforward to show that 

Z (E, A, d) = (27rm) A / 2 (7r/C) (A - 1)/2 • 2b A 1 / 2 ■ (\[ \\ /2 

\k=i 

. { P-^ + _g^ , e[g _^ +crf2)] 

_ y C H \ E - A(V + Cd 2 )- CR 2 ] A ~ m - 3/2 

^ n m\T(A-m- 12) ' J 

m=0 v ' ' 

■ Q[E — A(V + Cd 2 ) - CR 2 ] (8) 
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where R has been defined above. This expression shows explicitly that Z (E, A, d) 
has an algebraic dependence on E starting from a threshold energy E = A(V + 
Cd 2 ). There is no sign for an exponential behaviour which would indicate the 
existence of a plateau in the (T = P _1 ,E) representation. The fact that Z is a 
monotonic increasing algebraic function of E is confirmed by numerical tests of 
expression (8), see Fig. 6. This confirms and comforts the results obtained through 
simulations in section 2. 

The analytic results can easily be extended to a 3-dimensional system since, 
because of the choice of the two-body potential, the potential contribution factorises 
in the 3 space dimensions (x,y,z). Hence the expressions are formally similar 
to those obtained in the 1-dimensional case, except for the eigenvalue spectrum 
{Xk',k = 1, . . . ,A} which is affected by the dimensionality of the total space because, 
in the nearest-neighbour interaction approximation the coordinance goes over from 
2 in the 1-dimensional case to 6 in the 3-dimensional one. The partition function 
reads now explicitly 




Z (E, A, d) = (2nmf A / 2 ■ (2bA 1 / 2 ) 3 ■ (n/C)^/ 2 
■[[E- A(2V + 6Cd 2 )] 3A ~ 5/2 -Q[E- A(2V + 6Cd 2 )] /T(SA - 3/2) 

(A-3)/2 r mn2m „. r .„ 

_ 3 y , , , , r \E - A(2V + 6Cd 2 ) - CR 2 ] 3A "™" 5 / 2 

^ m\T(3A - 3/2 - m) ' J 



m=0 



• G [E - A(2V + QCd 2 ) - CR 2 ] 

(A-3)/2 

_|_ 3 \ 

^ mlm'\T(3A- 3/2 - m - m') 

m,m'=0 v ' ' 

. [ E - A (2V G + 6Cd 2 ) - 2C R 2 ] 3A - m - m '- 5/2 
■S[E- A(2V + 6Cd 2 ) - 2CR 2 ] 

(A-3)/2 ^ m+rn ' +rn " ff2(m+m'+m") 

^ mlm'lm"W(3A - 3/2 - m - m' - m") 

m,m,m"=0 
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. [ E - A(2V Q + 6Cd 2 ) - 3CR 2 Y A - m - m " m ~ b/Z 

■ 6 [E - A(2V + 6Cd 2 ) - 3CR 2 ] (9) 

Here R is given by the same expression as (5). 
The eigenvalues are 

n / 2m 27rj 2nk\ 

\,j,k = 2 I 3 - cos — - cos — - cos — I 

where (i, j, fc = 1, . . . , JV) and N = A 1 / 3 . 

The correlation between E and T is shown in figure 7. It is similar to the one ob- 
tained in the one-dimensional case. In practice, in both cases, the relation between 
E and T looks like that of an ideal classical gas T = 3/2E, but with a different 
slope because of the presence of an interaction. 

4. Discussion and conclusions 

We have worked out the properties of an excited system of nucleons with re- 
spect to its energy. We introduced a classical model which is aimed to describe a 
fragmenting nucleus considered as a disordered system. In this model the nucleons 
occupy finite phase space cells. We have shown in former work that this approach 
reproduces very successfully the fragment size distributions and related observables 
predicted by bond percolation models and experimentally observed, hence the prop- 
erties related to a second order phase transition corrected for finite size effects. 

Our interest in the thermodynamic properties of the system was triggered by 
recent experimental results which indicate that there may exist a first order phase 
transition in the fragmentation process. The coexistence of phase transtitions of 
first and second order could be explained by the fact that first order transitions 
become second order transitions at the critical point. The universality properties of 
the fragment size distributions may reflect the existence of events which lie at the 
critical point and in its neighbourhood. We used our model in order to construct 
the micro canonical partition function, taking into account the short range nuclear 
and long range Coulomb interactions between all the nucleons. It comes out of our 
calculations that we do not find the plateau of constant temperature in the caloric 
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curve which characterises a first order transition. This is the case for numerical 
simulations with different nuclear interactions and also for a simplified analytical 
approach which leads to a steady increase of the temperature with increasing energy. 
In numerical simulations we faced the problem that part of the generated config- 
urations do not correspond to thermal equilibrium, since we did not select them 
by means of an appropriate Metropolis Monte Carlo procedure. This is the reason 
for which we implemented an analytical approach which implies thermodynamical 
equilibrium by construction. In order to ensure an analytical expression for the 
partition function which shows explicitly the behaviour of the caloric curve, we in- 
troduced a simplified two-body potential. This potential reproduces the qualitative 
behaviour of the realistic interaction used in numerical simulations. 

The reason for the failure of the present model with respect to the presence 
of a phase transition are not clear. We suspect that the concept of cells in space 
introduces a rigid structure which leads to a system closer to a disordered crystal 
than to a liquid which may vaporize when the excitation energy is high enough. 
One may of course relax this constraint and possibly, at the same time, select those 
configurations which correspond to equilibrium. But then the new approach should 
also reproduce the fragment size distributions which are so nicely obtained in our 
present approach. 

It is our aim to investigate this point in the future and confront it with the 
outcome of cellular automata techniques which introduce a dynamical description 
through an explicit time-dependence. 
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Figure captions 



Fig. la: S(E) = in Z(E, A, d) as a function of E/A for A = 216 and d = 3 fm. 

Fig. lb: (dS/dE)^ 1 as a function of E/A (V = Ad 3 is the volume of the system). 

Fig. 2: (dS/dE) v as a function of E/A for different values of d (V = Ad 3 is the 
volume of the system) ; the dashed line corresponds to the case where S 
is averaged over different values of d (see text). 

Fig. 3: S(E) and (dS/dE)^ 1 as a function of E/A for d = 3 fm with the interaction 
taken from ref. [20]. 

Fig. 4: Evolution of the partition Z(E) and fragment multiplicities as a function 
of E/A. See text. Z(E) is represented in arbitrary units. 

Fig. 5: Comparison between the exact (4a) and approximate (6a) expressions of 
the canonical partition function. The quantities which are represented cor- 
respond to the product of error functions in (4a) and the term in brackets 
in (6a). See text. 

Fig. 6: s(E) = £nZ(E, A, d)/A and temperature T as a function of E/A for d = 
3fm and A = 125 (1 dimensional system). See text. The parameters are 
the same as in Fig. 5. 

Fig. 7: Same as Fig. 6 for a 3 dimensional system (A = 5 3 ). 
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